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Abstract 

The one-dimensional quantum harmonic oscillator problem is examined via the Laplace trans- 
form method. The stationary states are determined by requiring definite parity and good behaviour 
of the eigenfunction at the origin and at infinity. 
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I. INTRODUCTION 



The quantum harmonic oscillator is one of the most important systems in quantum me- 
chanics because it can be solved in closed form and its solution can be useful as approxima- 
tions or exact solutions of various problems. The harmonic oscillator is usually solved with 

n 

the help of the power series method [1|, by the algebraic method based on the algebra of oper- 
ators 12], or by the path integral approach [31. In recent times, the one-dimensional harmonic 
oscillator has also been approached by the Fourier transform method [4j-[i|. Another oper- 
ational method useful in quantum mechanics is the Laplace transform method. The Laplace 
transform method was used at the first years of quantum mechanics by Schrodinger into the 
discussion of radial eigenfunction of the hydrogen atom Q , and more than forty years later 
Englefield approached the Schrodinger equation with the Coulomb, oscillator, exponential 



and Yamaguchi potentials [9]. More than twenty years went by and the hydrogen atom was 



again examined with the Laplace transform method 10] . For some years now the 1/x 111 ]. 
Morse 12j], N- dimensional harmonic oscillator 13j, pseudoharmonic and Mie-type [14| . and 
Dirac delta [15] potentials have been solved for the Laplace transform. 

In Ref. {9}, Englefield found the spectrum of the three-dimensional harmonic oscillator by 
imposing that the radial eigenfunction vanishes at the origin and by using the closed-form 
solution for the Laplace transform. This paper approaches the one- dimensional Schrodinger 
equation for the harmonic oscillator with the Laplace transform method following the recipe 
proposed by Englefield [9f. Nevertheless, we do not use the closed-form solution for the 
Laplace transform. In addition, we enlarge the class of problems to include eigenfunctions 
satisfying homogeneous Neumann condition at the origin. The main features of our approach 
are as follows. After factorizing the behaviour at infinity, the second-order differential equa- 
tion for the eigenfunction ip (x) transmutes in a nonhomogeneous first-order equation for 
the Laplace transform The closed-form solution for $ (s) is not necessary. After in- 

specting the singularities of the differential equation for $ (s), we make a series expansion 
about the appropriate singular point, related to the proper behaviour of ip (%) near infinity. 
The behaviour of the eigenfunction near the origin, related to the behaviour of $ (s) near 
infinity, implies to truncate the series. The root of the indicial equation is the quantization 
condition. The recurrence relation for the coefficients of the truncated series for $ (s) is 
solved and the inverse Laplace transform is identified with the Hermite polynomials. This 



procedure should be helpful for students in a course on mathematical physics. 



II. THE LAPLACE TRANSFORM APPLIED TO THE HARMONIC OSCILLA- 
TOR 

Let us begin with a brief description of the Laplace transform and a few of its properties 



161 ] . The Laplace transform of a function f(t) is defined by 

POO 

F(s) = C{f(t)}= / dte- st f(t). (1) 
Jo 

If there is some constant treK such that 

|e _<rt / (<) | < M, (2) 

for sufficiently large t, the integral in Equation ([[]) will exist for Re s > a and / (t) is said 
to be of exponential order. The Laplace transform may fail to exist because of a sufficiently 
strong singularity in the function / (t) as t — > 0. In particular 

f t a 1 1 
£ <i ttt 7 f = — TT) a > -1. (3) 

The Laplace transform has the derivative properties 

n-l 

£ {/w ( t )} = s n c {/ (t)} - J] (o) 

fc=0 

(4) 

£{*"/(*)} = (-l)"F<»> (a), 

where the superscript (n) stands for the n-th derivative with respect to t for (£), and 
with respect to s for (s). If near a singular point so the Laplace transform behaves as 

s^sq — So) 



then 



fit) ~ Trfr^V ', (6) 



where V {y) is the gamma function. On the other hand, if near the origin / (t) behaves like 
t a , with a > —1, then F (s) behaves near infinity as 

F{s) ~ (7) 
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We are now prepared to address the quantum harmonic oscillator problem. The one- 
dimensional Schrodinger equation for the harmonic oscillator reads 

H 2 d 2 1 



+ -mw V — E \ ip(x) = (8) 



2m dx 2 2 

Because the harmonic oscillator potential is invariant under reflection through the origin 
(x — > —x), eigenfunctions with well-defined parities can be built. Thus, it suffices to concen- 
trate attention on the positive half-line and impose boundary conditions on ip at the origin 
and at infinity. Normalizability demands if) — > as x — > oo. Eigenfunctions and their first 
derivatives continuous on the whole line with well-defined parities can be constructed by 
taking symmetric and antisymmetric linear combinations of ip defined on the positive side 
of the x-axis. As x — > 0, the solution varies as x s , where 5 is or 1. The homogeneous 
Neumann condition ( dip/dx \ Q = 0), develops for 5 = but not for 5 = 1 whereas the 
homogeneous Dirichlet boundary condition (ip (0) = 0) develops for 5 = 1 but not for 5 = 0. 
The continuity of if> at the origin excludes the possibility of an odd-parity eigenfunction for 
5 = 0, and the continuity of dip/dx at the origin excludes the possibility of an even-parity 
eigenfunction for 5=1. Thus, 



ip (x) ~ x 5 , 5 

rr->0 



for if} even 

(9) 

1 for if) odd. 



On the other hand, the normalizable asymptotic form of the solution as x — > oo is given by 

iP(x) ~ e -™^7(2fi)_ (10) 

X—+OC 

This behaviour invites us to define £ = mux 2 /H in such a way that Equation is written 



as 

n 



*4+~+M^(0=o, (ii) 



' di 2 2d£ 4 
where k = 2E / (fvjS). The solution for all £ can now be written as 



^(0 = (0 e~ f/2 , (12) 



where the unknown 0(£) is solution of the confluent hypergeometric equation 17] 



^ 2 0(O , (h ^#(0 ,,,, n n ^ 

e- 5 ^- + (6-O- 5r -a^(O=0, (13) 
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with a = (1 — k) /4 and b = 1/2. Because the asymptotic behaviour of if) (£) is given by 
exp(— £/2) as £ — ?• oo, one has to find a particular solution of f JT3|) in such a way that (£) 
tends to infinity no more rapidly than exp (a£^j , with (3 < 1 and arbitrary a, for sufficiently 
large £. This occurs because — £/2 — >■ — £/2 as £ — >• oo. This condition, added by the 
fact that (£) varies near the origin as ^ <5/ ' 2 , ensures the existence of the Laplace transform 
of 0(0- 

Denoting $(s) = £{0(£)} and using the derivative properties of the Laplace transform 
given by Eq. (j3J), the transform of Eq. (Fl3|) furnishes the nonhomogeneous first-order 
differential equation for $ (s): 



. iN /3 fc + 3\ . . 0(0) 



(14) 



Notice that s = and s = 1 are singular points of this differential equation. To make use of 
the property of the Laplace transform near a singular point, Eqs. and ([H]), and taking 
into account the asymptotic behaviour of (£) near infinity, we try a series expansion of 



$ (s) about s = 



18|: 



Referring to Eqs. © and 



(s) = t 
J), we find 



£4 

i=o 



>> ^ 0. 



(15) 



$( S ) 



s—>oo 



^+ for -0 even 



for V odd. 



(16) 



Thus, the series terminates at j — n in such a way that 

1 for if) even 



n 



(17) 



and 



n + 3/2 for -0 odd 
(0) for if) even 



(18) 



arbitrary for if) odd. 

Beyond to make 0(£) to behave as as £ — >■ oo, another important consequence of the 
term s""^"' in Eq. ffl5|) is to give rise to the quantization condition u = (k + 3)/4 flo| . 



Thus, 



2n + 1/2 for ip even 



E n = fku x < 



(19) 



2n + l + 1/2 for ip odd. 
Inserting Eq. (TT5"j) into Eq. (j!4j) we obtain the recurrence relation 

n — j + 1 /2 for ip even 



» 



» 



x < 



J + l 

Inspection and induction yields 



O) (n)( — 1) J 
c) ' = c v ..' X < 



(20) 



ri — j for odd. 



rS-j+i/2) for ^ even 



(21) 



L (n-i)! 



so that 



$n (S) = t#° X < 



for ^ odd, 



r (n + 1/2) E"=o jirCn-Z+i/^"-^ 1 for ^ even 



™ ! EJ=0 j|(nJ)l5"-j+3/2 fOT ^ ° dd - 



(22) 



j!(n-j)!s""-''+ 3 / 2 

Using Eq. ((2]) and inverting the Laplace transform term by term, we can reconstruct <p n (£): 



<Pn (0 = X 4 



r(n + l/2)E? =0 j -,r(n4;T)r( 



- — — tor ip even 



n-j+l/2) 



(23) 



n - ^'=0 j!(n-j)!r(n-j+3/2) 

Using Legendre's duplication formula [17| 



6 n (£) turns out to be 



for ip odd. 



r (z + 1) r (z + 1/2) = 2- 2z y^T {2z + 1) 



(24) 



(n) 

n (O = ^X < 
'7T 



r (n + 1/2) £™ =0 ( 71^^)1 for V> even 



(25) 



, v n (-1)^(2^) 2 " +1 - 2J 
U - 2-*>j=0 j\(2n+l-2jV. 

Then, using the formula for the Hermite polynomial 17 ] 

-I) 3 ' (2y) n ~ 2i 
j\(n-2j)\ 



for ^ odd. 



(26) 



j=0 
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where [n/2] denotes the largest integer < n/2, the eigenfunction is found to be 

, (B) , . [ FJ W0 lH ^(VW x ) for ^ even 



^(x) = - 7 =exp(-— ,)x 



(27) 



^ for ^ odd. 

Because if n (— y) = {—l) n H n (y), the solution can also be expressed as 



E„ = ( n+ ^ ) /)^. ;/ = 0. I . 2. . . . 



(28) 



where A, are normalization constants. 



III. CONCLUSION 



We have shown that the complete solution of the one-dimensional quantum harmonic 
oscillator can be approached via the Laplace transform method with simplicity and ele- 
gance, even if the eigenfunction does not vanish at the origin. Englefield's recipe concedes 
to explore asymptotic expansions of the eigenfunction and its Laplace transform to gain 
information from the Schrodinger equation even if it is not possible to solve it in closed 
form. The discussion presented here may apply to any other problem that, after factorizing 
the behaviour at the origin and at infinity, reduces to the confluent hypergeometric equation 
such as the pseudoharmonic, Coulomb and Kratzer potentials. 
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